Bakamjian-Thomas mass operator for the few-nucleon system from chiral dynamics 
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We present an exploratory study consisting in the formulation of a relativistic quantum mechanics 
to describe the few-nucleon system at low energy, starting from the quantum field theoretical chiral 
Lagrangian involving pions and nucleons. To this aim we construct a Bakamjian-Thomas mass 
operator and perform a truncation of the Fock space which respects at each stage the relativistic 
covariance. Such truncation is justified, at sufficiently low energy, in the framework of a system- 
atic chiral expansion. As an illustration we discuss the bound state observables and low-energy 
phaseshifts of the nucleon-nucleon and pion-nucleon scattering at the leading order of our scheme. 
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I. INTRODUCTION 



The quantitative understanding of the dynamics of few-nucleon systems, has reached such a high degree of accuracy 
that the need of including relativistic effects has been put forward Starting from a non- relativistic setting, for 

^ . a given set of particles, it is possible to put constraints on the mutual interactions in order to fulfill the requirements of 

■ relativity. This problem has been posed in full generality by Dirac in 1949 [3], and developed among others by Krajcik 
f--^ , and Foldy [1]. The relativistic corrections to the potential ("drift potentials") can take the form of multinucleon forces 
fS| ' and are currently being investigated quantitatively Q; in some cases they are being used to help in the resolution of 

[ persisting puzzles in the few-nucleon physics, like the neutron-deuteron Ay problem Q. 

■ Since the pioneering work of Weinberg and collaborators [1, the nuclear interaction potential can nowadays 
I be thought of as derived from a (relativistic) quantum field theory based on the chiral symmetry, which is the low- 
. energy effective theory of QCD (see Rcf. [lO| for a recent review). A systematic perturbative expansion in powers 

of p/Ah (Chiral Perturbation Theory, ChPT), the ratio of typical momenta divided by the typical hadronic scale 
Ah '-^l GeV, is formally possible, whose convergence properties have to be checked case by case in actual calculations. 
2 ' For power counting purposes a non-relativistic reduction is usually performed (heavy baryon expansion, HBChPT). 
(-H , This happens at two levels: in the calculation of Feynman amplitudes and in the definition of the effective potential, 
formulated by Weinberg in the framework of old-fashioned (time-ordered) perturbation theory. While the first step 
is not strictly necessary (a relativistic covariant scheme for the calculation of the two-pion exchange NN potential 
exists [ll| which relies on the covariant version of baryon ChPT formulated by Becher and Leutwyler the second 
?H ' step seems unavoidable. Thus relativistic and chiral corrections get intertwined, despite the different character of the 
corresponding symmetries: Lorentz invariance is an exact symmetry of Nature, whereas chiral symmetry, although 
useful as organizing principle, is anyhow approximate. One might therefore want to treat relativity exactly. The 
question then arises whether it is possible to set up a framework in which relativity is built in from the beginning, 
and corrections come only from dynamics (in our case from higher chiral orders). To the best of our knowledge this is 
not done in current ChPT approaches to few-nucleon systems. By separating the two expansions, one would also be 
able to compare relativistic and chiral corrections, and check quantitatively the effectiveness of the commonly adopted 
combined chiral and non-relativistic expansion. Indeed such a test is only possible, in our opinion, by performing a 
complete relativistic calculation, within a given framework, and comparing to its non-relativistic limit. Moreover a 
fully relativistic setting would allow to describe particle production, contrary to non-relativistic quantum-mechanical 
treatments. 

The approach we take in the present paper is therefore to describe the few-nucleon systems with a relativistic 
coupled-channel wave equation, with simple transformation properties under the Lorentz symmetry, so that the 
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behaviour of e.g. a, N — N subsystem in relative motion with respect to other nucleons, as demanded by relativity, 
is easily accounted for. In this perspective we consider the proposal of Ref. [l^ (see also Ref. [3] for a similar 
point of view) consisting in a Bakamjian-Thomas construction formulated in the point-form of relativistic dynamics. 
This proposal has already been implemented for bound state problems in the framework of chiral quark models p^ . 
Compared to the above references we focus on nucleons instead of quarks and adopt a systematic ChPT-inspired point 
of view, instead of relying on specific models. We apply the above mentioned framework to nucleons interacting with 
pions according to the dictates of chiral symmetry, using the vertices of the chiral Lagrangian at the lowest order, 
and restrict ourselves to the 1 nucleon sector {it — N scattering) and to the 2 nuclcon sector (both the bound and 
scattering states of iV — N). By adjusting the four low-energy constants which appear at the leading order, we obtain 
a reasonable unified description of these systems at low energy. However, the emphasis in the present paper is not on 
the accuracy of the description (limited to a leading order treatment), but rather on the presentation of a relativistic 
formalism which can naturally extend beyond the threshold of pion production. 

The plan of the paper is as follows: a brief review of the adopted framework is given in the next section, in which 
the mass operator acting on the truncated Fock space is constructed from a Lagrangian density; in Section IIIII we 
give a power-counting argument to justify the truncation of the Fock space in the low-energy expansion in the case in 
which the interactions are restricted by chiral symmetry. Section HVl contains the application to the 1-nucleon sector, 
which yields the nucleon mass rcnormalization and a Lippman-Schwinger type equation for the n — N scattering. We 
adjust the low-energy constant appearing at this order to reproduce the S wave scattering lengths; the pion-nucleon 
axial coupling 5^ is fixed by the peripheral NN phaseshifts, obtaining a reasonable description of the phaseshifts at 
small laboratory momentum (less than 50 MeV). In Section |V] we consider the 2-nucleon sector. Two more vertices 
(contact interactions) arise at the lowest order and we adjust the corresponding coupling constants to reproduce the 
S-wave scattering lengths. The agreement with experimental phaseshifts is again reasonable, up to center-of-mass 
kinetic energies of 100 MeV. 

II. MASS OPERATOR FROM A LAGRANGIAN DENSITY 

A. Bakamjian-Thomas construction in the point-form 

The requirements of relativity are given by the Poincare commutation relations among the generators of the group, 
written in terms of the particle coordinates. This translates into specific constraints for the possible interactions to 
include in the generators. Dirac classified three different possibilities each one associated with a particular spacelike 
hypersurface left invariant by a subgroup of the Poincare group: for the instant-form the hypersurface is the hyperplane 
t =const., for the point-form it is the hyperboloid — x.^ = t^, for the light-front it is the hyperplane t + z^O. The 
generators associated with these hypersurfaces are said to be "kinematical" , and do not contain interactions. In the 
case of the point-form the Lorentz transformations arc kinematical and are the same as in the free case. Only the 
four momentum contains interactions, and the requirements of relativity arc simply written covariantly as 

[P^ P''] = 0, UaP^U^' = {A-Y^ P", (1) 

where the generators of the Lorentz transformation A are unaffected by interactions. The Bakamjian-Thomas con- 
struction [l6l.[r^ in the point-form consists in the definition, starting from the non-interacting Poincare generators, of 
a mass operator Mq = \/ PqPq^ and a four- velocity operator such that Pq — MqV^; one then adds the interactions 
only to the mass operator M = Mq + Mj, and reconstructs the interacting four- momentum as = MV^. Poincare 
commutation relations are then satisfied provided the interacting mass operator is a Lorentz scalar which commutes 
with the four- velocity V^. It is therefore particularly convenient to consider the "velocity states" [l^li^: these are 
linear combinations of multiparticle momentum states which are eigenstates of the four- velocity operator. They have 
the nice property that all the particles transform with the same Wigner rotation under a Lorentz transformation. 
Starting with (non interacting) n-particle states \pi(yi,P2<^2, ■■■iPnO'n) with individual four-momentum pi and spin 
projection at, one defines the internal momenta by going, through a canonical boost Bc{v), to the ccnter-of-mass rest 
frame, 

k, = B:\v)p,. (2) 

By definition ^ k; = 0. Bc{v) is a rotationless boost which transforms the system from its rest frame to total velocity 
V = Pi / , where Wk is the relativistic energy of a free particle with three-momentum k. A velocity state is 

obtained from a multiparticle momentum state defined in its rest frame after a boost to overall velocity v by means 
of B,[v), 



ju, ki, cTi, k„, (T„) =Ub^ (v) 1^1 : ^1 ■>•■•■) ^n; ^n) ■ 



(3) 
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With this definition, using the Lorentz invariant normahzation for the momentum states, the velocity states are 
normaUzed as follows, 

(^,ki,ai,...,k„,a„|z;',k^,a;,...,k-„,a;) ^(27r)3" ^^y„ ^-""^ voS^iv - V) ^ - k' (4) 

A Bakamjian-Thomas construction in the point-form is thus accomplished in practice by defining the mass operator 
to be diagonal in the four velocity , which is conveniently expressed in terms of velocity states as 

(i;',k>^|Af,|«,k„cT.) = {v\ k>^|H(0)b, k,,aM27r)H\V - v) ^^I^^ , (5) 

m and m' being the initial and final relativistic energies, m = X^jWk;, m' = X^i'^k'. • In this equation 7i(x) is the 
density of the interaction Hamiltonian, which is a Lorentz scalar. Therefore, the above definition is frame- independent, 
in view of Eq. Q , and in principle any rotationally- invariant combination of the momenta ki and k- would be allowed 
to appear. H{x) will be taken as a sum of vertices constructed as local products of field operators. The structure 
function / is introduced in order to compensate for the neglect of the off-diagonal terms in the four velocity, and to 
regulate the integrals as well. From its definition the structure function / is dimensionless. From now on we will take 
for / a real symmetric function of its arguments, further specified as a Gaussian function centered around zero with 
cutoff A, times an additional cutoff function ^ of the relativistic invariants, which may be needed in order to regulate 
the integrals. 



/(m, m') = exp 



r\2 



m') 



2A2 



t (6) 



The cutoff A, within the effective theory implementation of the present approach, should be understood as the short 
distance scale at which new degrees of freedom start to become relevant. One such structure function is understood 
for each vertex of the interaction Hamiltonian, and the requirement that physics at low-energy is independent of the 
cutoff (provided it is large enough) should fix the running of the coupling constants with A. Compared to Ref. 
we have introduced in Eq. ([5]) a different normalization for the matrix elements of the interacting mass operator, in 
order to properly match, in the case when v ~ v' and m ~ m' , with the quantum field theoretical result [l9| 

= / ^'^^^(^(^) - (7) 

where in the point- form F{x) = and the factor rr? appears as a Jacobian in passing from the overall momentum 
conserving (5 function to the velocity conserving (5 function. 

B. A simple toy- model 

As an illustration of the general setting, we start by examining the simple example of a scalar nucleon field \E' 
interacting with a pion field </), where the interactions are provided by a Hamiltonian density of the form 7i(a;) = 
f/4'^^(a;)5'(a;)0(a;). (Here and in the following all products of field operators are understood as normally ordered.) 
Creation of nucleon-antinucleon pairs is neglected and a truncation of the Fock space to a given maximum number of 
pions is considered from the beginnning. In the 1-nucleon sector, truncating the states containing two or more pions, 
the mass operator takes the form 



+1 



where m^v is the physical nucleon mass, and is the relativistic 1-nucleon + 1-pion free particle energy. The 

counterterm 55;°" is needed for the mass renormalization. Due to the form of Ti{x), the interactions show up as off- 
diagonal entries in the mass operator. The nucleon mass renormalization and pion- nucleon scattering are described 
as eigenvalue-eigenvector problems for this mass operator. For instance, for the eigenvalue m^r, the physical nucleon 
mass, one finds an equation for the counterterm 

6l-^=g^K^{D,+,-mN)-'K, (9) 
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FIG. 1: Diagram contributing to the nucleon mass renormalization. 

with Di+i = LUk + uj^, having defined 



uji,= yjm% + k^, u;Z = VW+^. (10) 
and the operator K, which connects 2-particle to 1-particle states, 

(i/, k, -k\9K\v, 0) ^ / i^'^ k, -k|vI/t(O)vI/(O)0(O)|i;, 0)vo{27r)^6^{v - v') 



z;o(27r)-^^^(v-v'). (11) 



The superscript refers to the sector of the Fock space with baryon number 1. The mass operator commutes with 
the baryon number, and there is the freedom to choose a different structure function / for each sector of the Fock 
space. Taking the expectation value of the above equation between 1-nucleon states and inserting a complete set of 
velocity states in the subspace of 1-nucleon + 1-pion states one arrives at the nucleon mass renormalization due to 
the "pion cloud" , 



= ^ i^ r ^"-'7--'-^' . (12) 



f rf^k 1 |/W(m^,o.k + ^£) P 
2mjv J (27r)3 4:UJkLu^ uj^ + lu^ - tun 

Eq. (fT2|) determines the counterterm Sl^'^ for each choice of the cutoff A and coupling constant g, and corresponds 
diagramatically to the process shown in Fig. [TJ 

In the 2-nucleon sector, an analogous equation describes the dcuteron, 

{D2 + SD ^2 + g'-KHmo - D2+i)-'Kcf,^ - mocj)?, (13) 

where <j>2 is a state vector in the subspace of 2-nucleon states, and the operators D2 and D2+1 are respectively the 
relativistic 2-nucleon and 2-nucleon -I- 1-pion energy. As in the 1-nucleon sector, a counterterm 62^'^ is introduced in 
the corresponding diagonal element of the mass operator, in order to properly renormalize the 2-particle states. The 
kernel K connects 2-particlc to 3-particle states and has matrix elements between velocity states 



(«',ki,k2,-ki-k2|gA>,q,-q) = ,g^^^4=%^z;o(2^)3j3(^_^,) (^4) 
x(27r)3 [2c^ki<53(ki - q) + 2uo^,5^{ki + q) + 2uo^j''{k2 - q) + 2uo^J^{k2 + q)] , 



with m = 2cjq and m' = cjki +^^1:2 +'^ki+k2' ^'^'^ what concerns the covariance properties, the structure function 
/(^^ for the 2-nucleon sector can be chosen different from f^^\ However, as we will sec in subsection III CI in order to 
have a consistent renormalization we have to choose /'^^■' = /'^-'. By left-multiplying Eq. (fT3|) with the bra (w,k, — k| 
representing a 2-nucleon state with four velocity v and relative momentum (in the centcr-of-mass system) 2k, one 
arrives, after insertion of a complete set of states in the subspace of 2-nucleon + 1-pion states, to an eigenvalue 
wave equation for the center-of-mass wave function 0^(k) = {v = (1, 0), k, — k|0^). Using Bose symmetry [(fiJ^ik) = 
(/)|^(— k)], the bound state equation becomes 

(2t^k + <5r(k))<^f(k)+2c.kA(k)0f(k)+ y'^i?(k,q)0f(q)=mz,0f(k), (15) 



with 



d^q g 



2 



/(2)(2c^k,C^k+Wq + C^k 



A{k)= / ^ ^1^+q^-q^ (16) 
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FIG. 2: Diagram corresponding to the disconnected kernel of the two-particle scattering. 



and 



i3(k, q)=: , ., hq^-q. (17) 

The term proportional to A{k) represents a wave function renormalization of the two-nucleon state: it describes 
diagrams in which the nucleon hnes arc disconnected and dressed with pion loops, cfr. Fig. [21 

It is possible to resum such diagrams and obtain a Lippmann-Schwinger equation for a reduced (renormalized) 
amplitude containing only a connected kernel. Alternatively, one can choose the countcrterm Jj"" so as to cancel the 
disconnected kernel, ^""(k) = — 2a'k^(k). 

Correspondingly, the NN scattering is described by the Lippmann-Schwinger equation, 

■^2 = c^2+9'' Ws-D2Y^ [^/:s- D2+iY^ K *2, (18) 

conn 

where -^/s = i? + ie is the scattering energy, and (f>2 is an eigenstate of the free mass operator in the two-particle 
subspace. The above equation can be equivalently written in terms of the scattering amplitude, T(q, k) defined as 

T{ci,\^)=g^v,^,-q\KYVs-D2+iY' K \^2W, (19) 

conn 

where k denotes the incident three- momentum of the interacting state ^2- Inserting a complete set of velocity states, 
and omitting the three-delta over the velocities (such factors will always be present because they appear in every 
interaction vertex) the LS equation takes the form 



T(q,k) = T/(q,k)+ /^P(^ 

where the potential. 



d3p y(q,p)T(p,k) ^20) 



F(q,k)-.g2(z;,q,-q|/vt[yj_^2+i] 'k k, -k) = B(q, k), (21) 

conn 

consists only of the connected kernel B in Eq. (jl7[) . with the substitution mo ^/s. 

C. Consistency of the renormalization procedure 

The renormalization of the 2-nucleon lines describing NN scattering, realized by the choice of the counterterm 
(52'^"(k) = — 2a;k^(k), and of the 1-nucleon line, Eq. p2|) . correspond to the same physical processes, as can be seen 
by comparing Figs. [T] and [2l Physical considerations would require that, when the two nucleons are far apart and at 
rest, their energies should be renormalized as their respective masses. This implies the condition 

(5r'(0) = 25^°'^, (22) 

which can be regarded as the manifestation of the cluster decomposition principle in the simple case of two particles. 
The equation to fulfill is therefore 



(27r)^ [ Smj^ujqUj'^ tun — UJq — cj^ j 2771^ J {27t)^ 4u;qa;^ ujq + oj^ - mjsf 
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where we have replaced in Eq. ([T6)) by ^/s = 2m n, since we are considering the case of two widely separated 
nucleons at rest. It is seen that our choice of the structure function, J*-^^ = Z*-^^ = / depending on m — m' as in Eq. ([6|. 
independently of the baryon number sector, satisfies the requirement of a consistent renormalization procedure. Notice 
that this would not happen had we chosen the original formulation of Ref. : the crucial point was the inclusion of 
a different normalization for the matrix elements of the interacting mass operator, Eq. ([5]), which in turn was dictated 
by a proper matching to the quantum field theory. The cluster decomposition principle, satisfied by local quantum 
field theories, could in general be violated by a truncation of the full quantum field theory to a relativistic quantum 
mechanics. In view of the above consideration, we will from now on drop the superscripts and use the same structure 
function / for all sectors of the Fock space. 



III. EFFECTIVE THEORY IMPLEMENTATION 



Having described the general features of the construction of the interacting mass operator from a vertex Lagrangian, 
we now proceed to make full use of the constraints given by chiral symmetry. Most importantly, the Goldstone theorem 
requires that the coupling between pion and nucleons be of derivative type (suppressed at low energy), and this will 
in turn provide a power-counting justification for the truncation of the Fock space, since the creation of pions brings 
more and more powers of momentum. A lowest order interaction Lagrangian which respects chiral symmetry is 



(24) 



In this expression ip is to be understood as an isospin spinor and the pion field tt = t^tt" as the Goldstone boson 
SU(2) matrix, representing the coordinates of the coset space of chiral symmetry breaking. In addition to the 
pion-nucleon-nucleon vertex, at lowest order one has to consider other vertices such as the Weinberg- Tomozawa or 
the nucleon-nucleon contact interaction, which comes into two independent operators, but let us for the time being 
concentrate on the above vertex. It will contribute to the mass operator as off-diagonal matrix elements. For instance, 
in the two-nucleon sector. 



( D2 K 

K D2+1 K 

K D2+2 

V 



\ 


( ^ 




( \ 




02+1 




02+1 




02+2 




02+2 




I ■■ ) 




I ■■ 1 



(25) 



where we have defined for ease of notation D2 = D2 + 52°^^. Here K connects states containing n and n + 1 pions and 
is defined as in Eq. ([5|) using Ti.{x) = —Ct^m- The multi-particle states 02+n are expanded in velocity states which are 
normalized as in Eq. Counting three- momenta as small parameters of order 0{p), the chiral power of the phase 
space element d2+n<t', such that 



(^2+n0|02+n)(0: 



'2+n| 



in the two-nucleon sector is 



d2+„0'-O(p'"+'). 



(26) 



(27) 



The formal solution of the above eigenvalue problem gives an equation for the two-particle component of the state 
vector. 



CS-D2-K 



1 



^/s - D2+1 - K- 



-K 



-K 



= 0. 



\/s— r>2 + 2 

One can Taylor-expand the denominator in the continued fraction (j28[) . and compare the different operators 

1 2 



{,^fs-D2)-K- 



1 



rs-D2+ 



-K -K 



1 



2+1 



K- 



I 



^- D 



-KK..., 



2+2 



(28) 



(29) 



contracted between two-particle states. Since we are interested in kinematical configurations in which ^/s — -D2 ^ 0{p), 
the first term goes like 



(2|Vi-^2|2')^0(p-2); 



(30) 
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FIG. 3: Connected diagrams contributing to the ttN scattering at the leading order, corresponding to the axial coupling (upper 
panel) and to the Weinberg- Tomozawa vertex (lower panel). 

the second term, inserting a complete set of three-particle states, 

{2\K \ K\2') ^ 0(/), (31) 

where we have used the fact that {2\K\2 -I- 1) - 0(p"^), while the third term, since (2 + 1\K\2 + 2) - 0{p-'^), 

-\ 2 



{2\K 

The hierarchy is 



1 



D 



2+1 



K ^ \ KK\2') ^ Oip'). (32) 



0(p-2)+0(/) + 0(p2) + ... (33) 



Therefore the inclusion of 4-particle states (2-nucleons -I- 2-pions) ^2-1-2 yields a contribution suppressed by two orders 
in the chiral counting compared to the one of 3-particle states (2-nucleons -I- 1-pion). This observation is the basis to 
justify the truncation of the Fock space. It is possible to convince oneself that the same mechanism applies for the 
most general vertex respecting chiral symmetry. In this paper we stick to the first order of the low-energy expansion. 

IV. ONE NUCLEON SECTOR 

The vertices to be considered at the lowest order in the 1-nucleon sector for the study of TriV scattering come from 
the Lagrangian 



= -^V57^75a^7r^ + ^^7'^K,5^^]^. (34) 



These vertices are part of the leading chiral Lagrangian, expanded up to terms quadratic in the pion field. From the 
argument presented above, the neglected terms start to contribute at the next to leading order. Only the nNN vertex 
contributes to the nucleon mass rcnormalization, which reads. 



^ " AF^ J (27r)3 4wkwJJ Wk + TON uji, + uj^ - TUN ' 

The subscript in / denotes that the form factor is associated with the axial ttNN vertex of the interaction La- 
grangian. As in HBChPT, the mass renormalization counts formally as 0{p^). 

For the description of ■nN scattering, one needs to take into account the truncation u p to 2 pions, because the 
-kN states are connected through the leading vertex to pure N states and to tt-kN states • This corresponds to 
the fact that the scattering contains the direct nucleon pole and the crossed one, as shown in Fig. [31 where the 
Weinberg- Tomozawa vertex is also included. 

Consider then the pion-nucleon scattering in the center of mass frame, with {q, r, ({k, s,j}) the initial (final) 
nucleon momentum, spin and isospin indices, and {— q, a} ({— k, 6}) the initial (final) pion momentum and isospin 
index. The kernel is written as 

QVia, ^^ ^ 9a /A(t^q + CU^,mjv)/A(mA,,LJk + ^k) (t^T^ 
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FIG. 4: Disconnected diagrams contributing to the nN scattering at the leading order. 



+ 



xu(q, r)f-/5u{0, (7)u{0, cr)^'' 75w(k, s) 

g\ /A(2Wq,a;q+k+COq+t^J)/A(Wq+k+t^k+W^,2Wk) 

2^q+k y^K+C.-)3K+^£)3 - ^q+k - ^5 - 

xu(q, r)|(f''75u(-q - k, cr)?2(-q - k, (T)/'75M(k, s) 

JH^T^ k. q q^ i2(q,r)(r+rMk,.), (36) 



where we have denoted by and fc'^ the pion 4-momenta; the first two terms come fi'om the direct and the crossed 
nucleon pole, while the third one from the Weinberg- Tomozawa vertex, with an associated structure function fwT- 
The Dirac 4-spinors are normalized so that ""(k, s)u{'k, s) = + m^. The isospin structures are proportional to 
6ab, ("isoscalar") that we will denote with a "+" superscript, and to ie°-^'^T^, ( "isovector" ) that we will denote with a 
"-" superscript. Analogously, the spin structures are proportional to 5rs and to q x k • cr,,^. The kernel can thus be 
written in the following operatorial form, 

5(q,k) = 5"'' [g+((7,fc,cose') +iqx k- cr/i+(g,fc,cos6')] 

-|-ie"^V= [g~ {q, fc, cos 61) + iq x k • cr/i" {q, k, cos 9)] . (37) 

Notice that at the same order there is also a disconnected kernel, corresponding to the diagram shown in Fig. H) 
Analogously to what was done in the simple model of the previous section, the disconnected contribution to the 
scattering kernel can be absorbed by a local counterterm (5™\ in the corresponding diagonal entry of the mass 
operator. 

Thus the corresponding Lippmann-Schwinger equation reads, 

I l.\ hi I f ^T'P^dp bl^{q,p)tl^{p,k) 

te±{q,k) =bt,±{q,k)+ / -f—^ — ^T-T"' (3^) 

J {2ttY Vs-^p-^p+«e 

where 

bj,±{q,k)= f dx{g\q,k,x)Pe{x)+qkh^{q,k,x)[Pe±i{x)-xPe{x)]} , (39) 



the functions {y ~ g,h) correspond to ?/^/^ = + 2y^ and i/3/2 — y^ ~y^ ^ and the subscript ± corresponds to total 
angular momentum J = ^±l/2. At the first order of the low-energy expansion, we can neglect the contributions of 
The equation is then solved for the off-shell amplitude t^ij.{q, k) by discretizing the energy domain of integration, 
for different choices of the c.o.m. energy -y/s. Phaseshifts are then calculated by putting the particles on energy shell, 

^2.5',^iE) ^ 1 _ 2m^^tl^ip,p) = 1 + 2zp//±(p), (40) 

where p is the c.o.m. momentum, E = y/ mj^ + p'^ + \/ + p"^ . In the non-relativistic limit, realized as tojv — *■ oo, 
our results coincide with the ones obtained in the heavy baryon formulation of ChPT [20| . at the same chiral order. 
Although no complete fit procedure has been performed, we have chosen the Weinberg- Tomozawa coupling constant 
so as to reproduce the 7=1/2 5- wave, while gA is fixed from the peripheral NN phaseshifts (to be discussed later). 
This is repeated for three values of the cutoff A = 300 — 400 — 500 MeV. In Fig. [5] the lowest partial waves are 
shown as function of laboratory three- momentum, compared with the experimental analyses (squares) [2l|. The 
bands represent the uncertainties obtained by varying A between 300 and 500 MeV. However this is only a lower 
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FIG. 5: (Color online) S and P waves of ttN scattering as function of laboratory momentum. The (blue) bands are obtained 
by varying the cutoff A between 300 and 500 MeV. Also shown are the experimental analysis [2l|] (squares), our results in the 
non-relativistic limit for A = 400 MeV (dashed) and the three-level HBChPT result (dotted) . 



bound on the theoretical uncertainty: in particular one should expect a larger uncertainty to come from the neglected 
higher orders, especially from the contribution of the NLO couplings c^'s which are known to be anomalously large 
because of the proximity of the A (2^. The description of P-wave phaseshifts is rather poor in the channel of the A, 
and reflects the lack of important physics not included in our leading order calculation. The dashed lines in Fig. [5] 
represent our result taken in the non-relativistic limit {rriN — > oo) for A = 400 MeV, for which the fit procedure has 
been repeated, and show that relativistic corrections to this leading order calculation are sizeable in the case of ttN 
scattering. The dotted lines represent the tree-level HBChPT calculation, which should be of comparable accuracy as 
ours. One should emphasize that the latter calculation is perturbative while ours takes into account, through the LS 
equation, the loops due to rescattering. These loops in the (perturbative) HBChPT framework come always together 
with subleading vertices, which absorb the cut-off dependence. The renormalization of effective field theories at the 
non-perturbative level is an intensively studied problem, currently the subject of controversies among experts [23j . 
The importance of rescattering effects is especially apparent for the / = 3/2 S'-wave. 
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V. TWO-NUCLEON SECTOR 

In addition to Ct^n, for the two-nucleon sector one has to consider, at the leading order of the fow-energy expansion, 
also a contact interaction Lagrangian C^n containing 4 nucleon fields and no derivatives. We therefore consider in 
this section 'H{x) = —C^^n ~ ^nn a-nd discuss separately the onc-pion exchange and the contact interaction. 



A. One-pion exchange 



The nucleon- nucleon wave equation can be written as Eq. (fT3|) . The two-particle component of the wave function 
in momentum space will carry in this case spin (r,r') and isospin indices, and in the center-of-mass system is 
written as 



(2wk + '5r(k))</)2(k), 



-I- 



ss' 



rr' •> 



where 



</'2(k). 



(1,0); k, r,i;-k, /,z'|02). 



(41) 



(42) 



In analogy to Eq. (|15p , the Lippman-Schwinger equation will again contain disconnected contributions of the kind of 
the function A and connected contributions of the kind of the function B, 



^CJ(k,q) 



rr' Ai' 



(k,q) + (27r)^J'^(k - Ci)6rs5r's'5^J5^,yA{^<i) - (q <^ -q, s ^ s', j ^ j') 



(43) 



with 



R (V ct\ — 



/yl(2Wk,IiJq+iiJq_k+'^k)/,4(l^k+l^q_k:+'^q>2Wq) 



A(k) 



\/s-Wq-Wk-l^q_k 

')7''75"(-q, s'] 

l/A(2wk,Wq+cjJ_^-|-a;k)l^ 



xu(k, r)7''75u(q, s)u(-k, r')7''75u(-q, s')plpln. 



XM(k, r)7^75u(q, r)u(q, r')7'^75u(k, r')pZpl, 



(44) 



where is the form factor entering the definition of the interacting mass in Eq. ([5]) corresponding to the first 
operator in Eq. ([34]), and we have denoted by p'^ the 4-momentum of the intermediate pion, p'^ = (Wq^j^, q — k). 
The disconnected contribution can be absorbed by the choice of the countcrtcrm (52°"(k) = — 2a'kA(k). Counting the 
involved three momenta as small quantities of order 0{p), the pion mass M^r ~ 0{p) and the energy denominators 
^ 0{p) we have that ^ 0{p^). The remaining connected kernel B (matrix in spin and isospin) reads explicitly. 



5(q, = So(q, k) [(<T . pi)^(<T . pi)^ - {rr ■ p^Yrirr ■ p,Y,\ \2Si Sf, - SiSj, 
with the vectors pi and p2 defined by 

k^q — q^k 



{k- 



-k;s ^ s';j ^ j'} , (45) 



Pi 

P2 



[wk -I- TO at] [wk + ruN 

k , q ■ 



Wk + niN UJq + TOJV 



and the scalar kernel 
So(q,k) 
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[uj„ + mpf\ [wk + niNl 



/A(2tJk, Wq + tJq_k + ^k)/A(t^k + Wq_k + ^<i^ ^tjq) 



cjk - t^, 



(46) 



(47) 



q-k 



B. Contact interactions 



The most general chirally invariant leading order NN Lagrangian can be written, after using Fierz reordering, as 
the sum of two terms 



(48) 
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where the notation for the couphng constants has been chosen so as to conform with the usual ones in the non 
relativistic expansion. The two vertices are diagonal in the Fock space basis that we have chosen, 



D2 + 5. 
gK 



rcn 
2 



Do 



C2 



<P2 
02+1 



92 
t>2+l 



(49) 



where we have denoted by C2 and C2+1 the contribution of the new interaction vertices to the interacting mass operator 
in the two-particle and three-particle subspace. According to the recipe, they contain in addition an associated 
structure function /. The contribution of C2+1 to the wave equation is of higher chiral order, therefore in our leading 
order calculation, we can neglect it and only include the operator C2 in the analysis. With this understanding the 
scattering equation projected in the two-particle subspace becomes 



292 



C: 



2<P2 



'K^ [V^-D2+i] gK\ 



= Vs02 



(50) 



The new vertices yield connected contributions to the wave equation. Let us denote this connected kernel, analogous 
to B(k, q)f.^/'^,f/ in the one-pion exchange potential, with the letter C. Its operatorial expression reads 



C(k,q) 



Cs 



1 - 



(cjk + m7v)(wq + ttin) 
qk 



LOk + TO7v)(a;q -I- niN) 
(k X q) • (Ti(k X q) • 0-2 1 /s(2tJq, 2wk) 



(k X q) • (cTi + (T2) 



-Ct 



(wk + TOw)(cJq + mw) + 



(q.k)2 



1 - 



(tJk + mAr)(tjq + rriN) 
(Wk - Wq)(wk + ttin) + k^ 



(tJk + "^iv)(wq + TOjv) 
(k + q) • CTi (k + q) • (72 
kq 



(Tl • (T2 



(cJk + m7v)(t^q - 
- t^k)(t^q + niN) 



- rriN) 

L „2 J 



+1 



(wk + mjv)(t^q + rriN) 
.(q- k) (k X q) • (cri + (T2) 



q • CTiq • cr2 



k • CTj^k • (72 



(wk + TOAr)(wq -t- niN) 
- {q ^ -q; 1 ^ 2}, 



q2k2-(q.k)2 1 /T(2c^q,2c^k) 

(Wk + m7v)(t^q + TOJV) J 8 



(51) 



where for the associated structure functions fs and Jt we take the same expression as in Eq. ([B]), with the choice 
^ =exp[— (g"* -I- A:'*)/(2A'*)]. Such additional cutoff for the contact interactions is needed to regulate the LS equation. 
It is chosen, analogously to what is done e.g. in Ref. [2^, so that the induced modifications are of higher order in the 
chiral counting (actually more than needed for our leading order calculation). The wave equation now reads 



w, 



B(k,q)- +C(k,q)-i^ 



</'2(q)5 



[^/i-2c^k] 02 (k) r 



(52) 



In the chiral counting i?(k, q) ^ C(k, q) ^ 0(1). The wave equation can be recast in the form of a Lippmann- 
Schwinger equation, from which phaseshifts are found by means of standard numerical methods (23 |. 

As already mentioned in the previous section, we have fixed the coupling from the ^/g wave, and the contact 
terms Cs and Ct from the '^S'l and ^5*0 scattering lenghts, requiring that 



a(3S'i) = 5.4 fm, 0(150) = -24 fm. 



(53) 



We have repeated this procedure for three values of the cutoff A = 300, 400, 500 MeV. The results are shown 
in Figs. Elini The shaded bands represent the variations with A and can be considered as the intrinsic theoretical 
uncertainty of our calculation. Also plotted, as dashed lines, is the result corresponding to A 400 MeV in the 
non-relativistic limit, realized as m^r — > 00. The differences make it possible to quantify the size of the relativistic 
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FIG. 6: (Color online) Spin singlet A''A^ phaseshifts as function of laboratory kinetic energy. The (blue) bands show the effect 
of varying the cutoff A between 300 and 500 MeV. The squares represent the Nijmegen phaseshifts [2^ . Also shown in dashed 
are the results in the nonrelativistic limit (mjv — » oo), computed with A — 400 MeV. 

corrections included in our scheme: they are always smaller than the effect of varying A. therefore, at least for 
the leading order of the chiral expansion, relativistic corrections are smaller than the neglected chiral orders. The 
peripheral waves are well described by the one-pion exchange, but the agreement with data is poorer for the lowest 
waves. Clearly some important physics is missing in our leading order calculation (see e.g. the discussion of Ref. [26j). 
This is the case, for instance, of the ^D^ wave: it is well known that two-pion exchange is very important for the 
description of this wave, and we do not include this process, not even implicitly, via subleading low-energy constants. 



The deuteron 



With all the parameters fixed by the scattering, for the chosen values of the cutoff A, we have computed the 
momentum space bound state mass and wave function (/if (k) in the ^Si-^Di channel. The configuration-space wave 
functions uq and U2 are found by Fourier transformation (cfr. Ref. Some bound-state observables are shown in 
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FIG. 7: (Color online) Spin triplet A'' A'' phaseshifts in the uncoupled channels {L = J). The legend is as in Fig. [6l 



Table [J In Fig. [TO]wc show the configuration-space deuteron wave functions, in which as before the bands represent 
the variations with A between 300 and 500 MeV. 



VI. CONCLUDING REMARKS 

In this paper we have considered the proposal of Ref. pli| of a point-form formulation of relativistic quantum 
mechanics consisting in the construction of a Bakamjian-Thomas interacting mass operator for the few-nucleon system. 
Taking the most general chirally invariant Lagrangian describing pions and nucleons, we have implemented the above 
construction in the framework of an effective theory, where the needed truncation of the Pock space is justified by 
a systematic low-energy power counting. This in turn relies on the restrictions that the underlying chiral symmetry 
imposes on the interaction vertices. By introducing a different normalization for the structure functions /, associated 
with each coupling of the Lagrangian, compared to the original proposal of Ref. [l3| . we obtain the correct matching 
with the point-form quantized field theory and show that the renormalization procedure is consistent with the cluster 
decomposition principle. 
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FIG. 8: (Color online) Spin triplet A^A'' phaseshifts and mixing angles in the coupled channels {L — J ± 1). The legend is as 
in Fig. H 



We have performed a complete leading order analysis of the TriV and NN systems by solving numerically the 
eigenvalue equations originating by the diagonalization of the mass operator, and obtained a rather coherent picture 
of the low-energy phaseshifts. Of course, especially in the ttN system, important physics, such as the effect of the 
A resonance, is missed at leading order, and the examination of higher orders is mandatory for a better quantitative 
description. By comparison with the non-relativistic limit of our framework we have assessed the size of relativistic 
corrections included in our scheme, and shown that in the NN sector they are always smaller than the accuracy of 
our leading order calculation: NLO chiral corrections are larger than our "all-order" relativistic corrections. Still, 
it could well be, provided the chiral expansion is convergent enough, that NNLO chiral corrections be smaller than 
our all-order relativistic corrections. This can only be checked by actual calculation. Moreover, the importance of 
relativistic corrections can depend on the observables: in Rcf. [6| they are found to be surprisingly large for the 
neutron-deuteron Ay . 

As a possible development then the next order can be analyzed, by including the 2-pion states and taking into 
account the next-to-lcading order terms in C. By examining the off-diagonal matrix elements of the mass operator, 
our framework allows to incorporate in a natural way the pion production channel whose understanding from the 
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FIG. 9: (Color online) Spin triplet A^A'' phaseshifts and mixing angles in the coupled channels {L — J ± 1). The legend is as 
in Fig. H 



field-theoretical point of view is only recently starting to emerge [28|. Furthermore, once the Lagrangian is fixed in the 
two-nucleon sector, the sector could be investigated. The point-form formulation, in which the boost operators 
have only a kinematical character, would allow to easily account for such relativistic effects as the boosting of the 
wave function of an interacting 2-nucleon subsystem, although the requirements of the cluster decomposition principle 
might be less simple to fulfill in this case. 
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A = 300 MeV A 


= 400 MeV A 


= 500 MeV Empirical (from Ref. [27]) 


Binding energy Ba (MeV) 


2.01 


1.93 


1.99 


2.224575(9) 


Asymptotic D/S state 77 
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0.018 


0.020 


0.0256(4) 


Matter radius (fm) 


1.86 
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1.971(6) 


Quadrupole moment Qd (fm^] 


0.12 


0.16 


0.18 


0.2859(3) 


D-state probability Pd (%) 


0.5 


1.2 


1.9 





TABLE I: Bound state observables for the chosen values of the cutoff A. 
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FIG. 10: (Color online) Configuration-space wave function of the S-wave (upper curves) and D-wave (lower curves) components 
of the deuteron. The (blue) bands show the effect of varying the cutoff between 300 and 500 MeV. 
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